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Abstract 

A model is presented in this work for simulating endogenously the evolution of the marginal costs of production of 
energy carriers from non-renewable resources, their consumption, depletion pathways and timescales. Such marginal 
costs can be used to simulate the long term average price formation of energy commodities. Drawing on previous work 
where a global database of energy resource economic potentials was constructed, this work uses cost distributions of 
non- renewable resources in order to evaluate global flows of energy commodities. A mathematical framework is given to 
calculate endogenous flows of energy resources given an exogenous commodity price path. This framework can be used 
in reverse in order to calculate an exogenous marginal cost of production of energy carriers given an exogenous carrier 
demand. Using rigid price inelastic assumptions independent of the economy, these two approaches generate limiting 
scenarios that depict extreme use of natural resources. This is useful to characterise the current state and possible uses 
of remaining non-renewable resources such as fossil fuels and natural uranium. The theory is however designed for use 
within economic or technology models that allow technology substitutions. In this work, it is implemented in the global 
power sector model FTT:Power. Policy implications are given. 



Keywords: 
depletion 



Global energy resources, Climate change mitigation, Economic Potential, Renewable Energy, Fossil Fuel 



X 



1. Introduction 

1.1. Energy-Economic-Environmental interactions 

The use of large scale models for exploring Energy- 
Economic-Environmental (E3) interactions is crucial for 
devising policy aimed at addressing coupled economic and 
environmental problems and achieve related policy goals. 
This is due to the fact that in such complex and highly 
correlated systems, while conceptual difficulties arise in 
attempting to understand the systems-wide influence of 
individual policies and regulations, significant complica- 
tions arise in the potential mutual influence between sev- 
eral such policies (Barker et al., 2007). This includes for 
instance the strong interaction between government sup- 
port for novel transportation technology and power sec- 
tor or land use management, and their very uncertain ef- 
fect onto global emissions, which depend highly on their 
timing, technology diffusion timescales and energy con- 
version efficiencies (as examples of differences in estima- 
tions of potential emissions reductions for the transport 
sector, see Pasaoglu et al., 2011, Takeshita, 2011, 2012, 
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van Vliet et al., 2011, 2010). It has been widely recog- 
nised that large expansions in modelling capacity are re- 
quired in order to better predict the likely result of com- 
prehensive policy portfolios, which should include combi- 
nations between top-down economic models and bottom- 
up technology models (see for instance Grubb et al., 2002, 
Koehler et al., 2006a,b). While common economic mod- 
els can evaluate the global demand for energy, transport, 
materials, goods and services, they generally do not rep- 
resent with much detail the way in which their supply is 
produced and at which costs, from lack of technology reso- 
lution, or none altogether. This generates for instance sig- 
nificant uncertainty over production efficiency and green- 
house gas (GHG) emissions. Meanwhile, bottom-up tech- 
nology models generally take demand values (energy, ser- 
vices, goods, etc) as given and therefore, although they 
are able to generate prices and accurate efficiency val- 
ues and emissions factors, they do not capture the inter- 
action between prices and demand (for details on these 
aspects for several existing models, see Edenhofer et al., 
2010, 2006)). Coupling bottom-up and top-down models 
generates the most powerful method to capture systems- 
wide and economy- wide coupled interactions, which are 
currently strongly required for devising sensible climate 
change mitigation policy (Koehler et al., 2006a). 

Energy flows, originating from natural resources, are a 
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necessary component for all sectors of the world economy. 
Although the economic output of the energy sector ac- 
counts only for a small fraction of the world gross do- 
mestic product (GDP), 1 changes in the prices of energy 
carriers have pronounced consequences on the output of 
most other economic sectors (see for instance Jones et al., 
2004)). 2 Since the price of energy carriers are reflected in 
the prices of goods and services originating from energy 
intensive sectors, such changes can lead to increased infla- 
tion, decreases in economic output and reduced paces of 
economic development. Many attempts have been made 
to capture such interactions between energy, the econ- 
omy and the environment in computer models (see for 
instance Edenhofer et al., 2010, 2006, and the various 
models reviewed). While many models of E3 interac- 
tions do not incorporate explicit representations of nat- 
ural resource use and depletion, or the physical limits 
to available energy flows, very few feature endogenous 
exploitation costs of non-renewable resources and none 
of them feature a particular emphasis on uncertainty in 
the economic or technical potentials of natural resources. 3 
For this reason, in previous work we defined a theoret- 
ical framework and built an extensive database with a 
resolution of 190 countries for limiting and tracking the 
use of natural resources in models of global energy sys- 
tems (Mercure and Salas, 2012), which, although adapt- 
able to any energy systems modelling framework, was de- 
signed for use in the model Future Technology Trans- 
formations in the Power sector (FTT:Powcr) (Mercure, 
2012a). FTT:Power is based on a theoretical framework 
for technology diffusion (Mercure, 2012b), and was in- 
tended to be integrated as a bottom-up component to 
the Encrgy-Economy-Environment model at the Global 
level (E3MG, for descriptions sec Barker et al., 2012, 2006, 
Barker and Scrieciu, 2010, Koehler et al., 2006a). 

Modelling energy systems realistically requires the rep- 
resentation of many complex interactions between differ- 
ent types of systems, which must respond to the economic 
climate and natural environment at every point in time. 
This involves modelling the behaviour of actors who influ- 
ence the working and composition of the technological mix 
within the system. Once this mix is defined, the require- 
ments in terms of energy resources are straightforward to 
evaluate. Global energy demand is strongly influenced by 
the price of energy carriers, 4 generating a feedback inter- 
action between the economy and the global energy system 
through demand and prices. Meanwhile, the cost of en- 
ergy resources influences the choice of investors in energy 
systems and thus the technology composition, as well as 
the cost of producing energy carriers. Therefore, a sec- 
ond strong feedback interaction exists between the global 
energy system and the natural environment through the 
exploitation of resources through demand and prices. As 
described earlier by one of us (Mercure, 2012a,b), the evo- 
lution of technology in most sectors, including the power 
sector, is well described by a coupled family of logistic 
differential equations that simulates transitions between 



energy technology systems, changes that arc driven by the 
trend of investor decisions, an approach supported by an 
extensive empirical literature (see for instance Bass, 1969, 
Bhargava, 1989, Grubler, 1990, 2012, Grubler et al., 1999, 
Marchetti and Nakicenovic, 1978, Morris and Pratt, 2003, 
Sharif and Kabir, 1976, Wilson, 2009). Meanwhile, the 
cost of producing energy carriers is influenced by that of 
natural resources, as well as and through components such 
as investment, maintenance, capacity factors and taxes or 
carbon costs, all of which should be considered when cal- 
culating the cost of electricity production, for which for 
instance the Lcvcliscd Cost of Electricity (the LCOE, see 
for instance ?), in the case of the power sector, is a good 
representation of the way investors evaluate technology 
costs (and in a similar construction for other sectors of 
technology). As argued in our previous work (Mercure, 
2012a, Mercure and Salas, 2012), the limitation of the ex- 
pansion of certain types of energy systems is well repro- 
duced by cost-supply curves, which track the increasing of 
the marginal cost of production of energy with increasing 
demand, through its influence into certain components of 
the the LCOE (e.g. fuel costs, capacity factors, investment 
costs, etc). 

Modelling energy flows from renewables and non- 
renewable resources entails large conceptual differences. 
Cost-supply curves have been generated for different 
types of renewable resources in works by de Vries et al. 
(2007), Hoogwijk (2004), Hoogwijk et al. (2004, 2005), 
van Vuuren et al. (2009), using the cumulative sum of cost 
rankings of the global number of potential sites for energy 
production by type (wind, solar and biomass energy) . This 
involves the assumption that these renewable resources are 
chosen and exploited in order of cost, starting with the 
most profitable development ventures. The cost-quantity 
availability of non-renewable resources such as oil and gas 
can also be expressed using a cost-quantity curve (as in 
Mercure and Salas, 2012, Rogner, 1997), which expresses 
a quantity available at a certain exploitation cost rather 
than a flow. Such a curve, however, is much more difficult 
to interpret in order to derive marginal costs, since tak- 
ing the assumption that consumption progresses in perfect 
order of exploitation cost is not reliable, and the gradual 
depiction of fixed amounts of resource means that the cost- 
quantity curve changes with time. In contrast, as apparent 
in the oil industry for instance, the exploitation costs of 
existing projects cover a wide range rather than a single 
competitive value, depending on the nature and quality 
of resource occurrences. This range is determined by the 
price of oil. There is thus a connection between the supply 
and the price of energy commodities, where higher prices 
enable production at higher costs, and therefore the ac- 
cession of larger amounts of resource at such costs. Mean- 
while, the demand for energy commodities may justify in- 
creases of prices, in order to enable production at higher 
costs, such that the demand is met by the supply, using 
ever more difficult and expensive methods, locations and 
types of resources (ultra-deep offshore drilling, arctic sites, 
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tar sands, oil shales, etc). However, high prices, as for in- 
stance generated by depletion and scarcity, may also be 
avoided by simply phasing out the use of certain types of 
high price commodities, replacing them by other types. 
Such substitutions actually stem from technology substi- 
tutions, which can become economical in the event of the 
price of some commodities increasing (e.g. replacing oil 
by coal for producing electricity, which entails phasing out 
existing technology, and therefore cannot be performed 
instantaneously 5 ). Thus, technology substitution has a 
strong long term influence as well on energy commodity 
prices. And conversely, the price of energy commodities 
has a strong long term influence on the technology mix. 

1.2. Oil price models 

Many models aiming at describing the dynamics of 
the price of oil exist in the literature (for instance 
Carlson ct al., 2007, Gallo et al, 2010, Michl and Foley, 
2007, Rehrl and Fricdrich, 2006, Reynolds, 1999). Models 
are mostly classified into two groups, based on whether 
they use the Hotelling Principle or the Hubbert Peak ap- 
proach (Reynolds and Baek, 2012). The Hotelling Prin- 
ciple generates the optimal extraction rate of a known 
non-renewable resource, where its price follows a rate of 
increase equal to either a social discount rate or an in- 
terest rate (Hotelling, 1931, Pcrman et al., 2003). Mean- 
while, the Hubbert Peak theory is an empirical observation 
that Hubbert made on historical US oil production data 
from 1900 to 1962 using a logistic trend (Hubbert, 1956, 
1962), for which a theoretical derivation was later devel- 
oped, based on probabilistic arguments concerning the rate 
of oil field discoveries, that generates the logistic (or more 
general) mathematical form of the observed trend. In this 
theory, while random drilling generates a rate of discovery 
which is proportional to the amount of undiscovered oil 
left in a geographical area, oil can statistically more easily 
be found near existing fields, generating a quadratic term 
in the probability function (the 'information effect'), lead- 
ing to a logistic differential equation (Rehrl and Fricdrich, 
2006, Reynolds, 1999). Both theoretical approaches can 
be criticised, for different reasons described below. 

While the Hotelling Principle does provide the optimal 
rate of resource extraction given a certain resource base, 
it is very unlikely that real resource extraction activities 
follow a path anywhere near optimal, for the following rea- 
sons. As Norgaard (1990) argues (using the 'Mayflower 
problem'), mineral extraction firms would need to know 
the total amount, location and quality of resources in the 
ground over which they have contracts, and would need 
to have perfect foresight, both required to enable them 
them to figure out how to follow the optimal pathway pre- 
scribed by Hotelling (see also Reynolds and Baek, 2012). 
There is no clear empirical evidence in the literature to 
justify the assumption that energy or mining firms actu- 
ally follow such patterns (see for instance the textbook 
Perman et al, 2003, pp. 527-532). 



Meanwhile, the extended Hubbert Peak theory gener- 
ates a supply that is entirely independent of demand, 
based solely on the rate of discoveries. Independent and 
individual Hubbert peaks have been observed for differ- 
ent types of oil occurrences (conventional US oil fields, 
deepwater, alaska, as described in Rehrl and Friedrich, 
2006). However, it has also been shown that OPEC be- 
haviour does not follow a logistic trend and has a reserve 
to production ratio of about 80 years, higher than the 
rest of the world (see below in section 2.3, as well as in 
Rehrl and Fricdrich, 2006). Thus it does not apply to all 
situations, and if it did, it would mean that oil consump- 
tion would follow a sum of rigid individual logistic func- 
tions, independent of oil demand, or alternatively with the 
demand and supply independent of the oil price. 

In contrast to this, for example, as the price of oil passed 
the economic threshold of 85-95$/boe (NEB, 2011) in early 
2007, intense activity started in the Canadian tar sands, 
which subsequently came to a standstill slightly later in 
2008 when the price dipped below that value again, and 
started again later in 2009 when the price increased again 
above a similar threshold. 6 This indicates how some ex- 
ploitation activities occur very near the competitive mar- 
gin, the production of which can be turned on or off, fol- 
lowing changes in demand rather than a Hubbert peak. 
Thus, the difference between a price dependent supply or 
demand and a sum of rigid Hubbert cycles is buffered by 
high cost resources at the edge of profitability and/or by 
OPEC monopolistic behaviour. Sorrell et al. (2010) states 
that 'Most of the worlds conventional oil was discovered 
between 1946 and 1980 and since that time annual pro- 
duction has exceeded annual discoveries', indicating that 
a significant fraction of the original oil resources in place 
have already been discovered. These are not necessarily 
under exploitation, and therefore discoveries and changes 
in the knowledge about the location of oil fields is appar- 
ently not the single determinant of oil demand and sup- 
ply. Finally, Hubbert peak theory ignores the dynamics of 
technological change and technology substitution, which 
influences the demand by transiting away from expensive 
resources. We thus argue in this paper that neither of 
these two strict approaches arc appropriate to project en- 
ergy carrier prices, and propose an alternative model that 
does not assume knowledge about resources, perfect fore- 
sight, or that is based strictly on rates of discoveries. 

1.3. New approach to modelling carrier flows and prices 

In this paper, we present a detailed description of a 
theoretical framework to treat the exploitation of non- 
renewable resources (stock resources henceforth) and the 
prices of their associated energy commodities. This model 
is designed to be used in conjunction with a database for 
natural resources as well as with a model of technology 
substitution. It is a general model that can be used to rep- 
resent any type of natural resource consuming systems, but 
is applied here within our power sector model FTT:Power. 
This model does not, however, treat the effect on prices of 
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speculation and hoarding, or of supply problems related to 
geopolitical events, making it unable to forecast short term 
price fluctuations that commonly occur in global markets. 7 
It does, however, project inevitable long term base price 
values below which the production sector cannot supply 
the demand, and therefore provides a lower bound for price 
values, equal to the marginal cost of production. It may 
additionally be argued that hoarding and speculation over 
futures can only ultimately lead to bubbles and cyclic fluc- 
tuations, since (1) speculation and hoarding cannot occur 
without storage space, and (2) storage space is finite and 
stored commodities must eventually return to the market 
for this activity to be profitable (i.e. storage space cannot 
expand indefinitely). Thus the artificial demand (i.e. de- 
mand unrelated to immediate consumption) generated by 
speculation is cyclic and evens out to zero in the long run, 
generating price fluctuations that oscillate. Price fluctua- 
tions observed occur at faster frequencies than the fastest 
possible rates of technological change, and are therefore 
seen as volatility (i.e. quantifiable risk). Therefore it is 
mostly long term price changes that are truly relevant in 
energy systems modelling, along with the cost associated 
to volatility, and thus the omission of hoarding and spec- 
ulation in this theory is not expected to affect the results 
significantly, and the cost of volatility can be included into 
the sum of fixed costs, added to the marginal cost of pro- 
duction. Similarly, only long term prices should be con- 
sidered in energy and climate policy. 

The results of this study provides insight into en- 
ergy policy-making, particularly in the context of climate 
change policy. Through the use of the concept of limiting 
assumptions, limits to the use of stock resources, as well as 
the limiting behaviour of commodity prices are explored, 
providing information on the range of possible futures for 
the use of stock resources, their prices and the potential 
scale of global emissions. These scenarios are price inelas- 
tic, however, and more realistic scenarios are given, ob- 
tained when connecting this theory with a model of tech- 
nological change. While this insight can be used directly 
in better understanding the current state of global stock 
energy resources, this model reveals its real value through 
its integration into FTT:Power, and in future work into 
the full FTT family of technology models interacting with 
E3MG. 

This paper progresses as follows. We first present a for- 
mal definition of a dynamic model for tracking the use and 
depletion of stock resources based on cost distributions. 
This model can be used with two types of assumptions, 
either an endogenous commodity supply given a commod- 
ity price (the forwards problem), or an endogenous com- 
modity price given an endogenous demand (the reverse 
problem). As modelling exercises, the model is used with 
our natural resource database (Mercure and Salas, 2012) 
in these two limiting assumptions cases to determine, with 
confidence ranges related to uncertainty in resource assess- 
ments, limits to potential future supplies or prices in situa- 
tions of fixed technology. The results of these two exercises 



are instructive as they provide a clear picture of the limits 
to stock resources. However, they do not generate realis- 
tic scenarios of resource use, since they use assumptions 
where prices are independent of the economy and of one 
another, and demand values are independent of prices. In 
order to create fully dynamic scenarios that include possi- 
ble technology substitutions, this model is integrated into 
FTT:Power. Two commodity supply and price scenarios 
are given where technology substitutions in the power sec- 
tor enable the energy system to avoid price escalations. 
The impacts of these considerations onto climate policy 
are discussed. 

2. Theoretical framework 

2.1. Flows from resources to reserves and the consumption 
of reserves 

A fixed energy carrier demand, when met by renewable 
energy resources, results in fixed levels of use of these re- 
sources. However, when met by stock resources, it results 
in particular rates for their depiction. The gradual de- 
pletion of stock resources generates gradual increases in 
their cost of exploitation, an effect which is due to the 
natural average tendency towards the exploitation of re- 
source occurrences with lower extraction costs first, and 
costs increase following their depletion. In our previ- 
ous work Mercure and Salas (2012), we calculated cost- 
quantity curves for fossil and nuclear resources. When 
allocating energy demand between all potential energy 
sources, dynamic rates of exploitation of stock resources 
emerge, which result in particular lengths of time for these 
resources to reach depiction, and which depend on the 
price of the associated energy commodities. We introduce 
in this section a simple mathematical model that generates 
a relationship between rates of exploitation of a given fi- 
nite resource base and the price of the associated commod- 
ity. This relationship is not functional however, but as we 
show, involves hysteresis and irreversibility, and therefore 
strong path dependence. The supplementary information 
provides additional mathematical details in order to un- 
derstand all of its properties, which are not shown here 
in order not to lose the reader into strongly theoretical 
considerations. 

In a hypothetical world with a perfectly efficient and 
competitive global energy market, fast rates of resource ex- 
traction and no speculation over the future value of these 
commodities, or monopolistic or cartel price markups, only 
the resources with lowest extraction costs would be traded 
at any time. In such a situation, resource consumption 
would proceed in perfect order of increasing extraction 
cost (or level of technical difficulty of recovery) and fol- 
low closely a cost-quantity curve. This is of course not 
the case, and in reality, owners of low cost reserves delay 
their extraction, forcing the exploitation of more expensive 
resources while low cost resources are not yet depleted. 
Thus, the demand for energy commodities is met by re- 
source types with extraction costs within a certain range, 
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Figure 1: Process of stock resource extraction, a) Original (current) resource distribution as a function of extraction cost (black curve) and 
probability function for resource extraction given the price of the resource commodity that delimit reserves out of resources (red curve), its 
maximum being equal to 1. b) Dynamic process of resource depiction as the price of the resource commodity increases. The black curves 
corresponds to the distribution of resources left after increasing amounts of time have passed, associated with increasing prices for the energy 
commodity, c) Flow curve and uncertainty range for the amount of resource unlocked by an increasing exogenous price of the commodity 
which allows the resource exploitation to proceed up to an upper limit marginal cost (function of time). The area between the red curves 
indicates the 95% confidence level region, while the blue curve represents the most probable curve, d Marginal cost of production and its 
uncertainty range of the commodity for an exogenous demand (function of time). The values diverge at depiction, which, due to uncertainty, 
can occur at different values of cumulative production. 



delimited by the price of energy carriers, which determines 
what is economic and what is not, giving a corresponding 
range of margins of profit for different resource types. For 
example, the price of oil determines which of known oil 
fields are deemed economic to exploit, and the remaining 
fields are reserved for a future in which a higher price of 
oil is expected. Additionally, however, the extraction of 
existing economic reserves is done over a length of time, 
and some reserves are kept untapped for a future where 
higher prices are expected. The upper limits for the ex- 
traction cost values that are still considered economical 
given the prices of energy carriers, in other words the cost 
of the most expensive resource exploited (i.e. the marginal 
cost), are defined by the differences between the prices of 
the commodities and the sum of all fixed costs such as 
those associated to transformation processes and trans- 
port, and the minimum profit margins that industries will 
consider. Increases in the prices of energy carriers enable 
wider ranges of natural resource types to be exploited, for 
instance low grade or unconventional fossil fuels, which are 
not profitable under low price conditions. 

This behaviour can be summarised by stating that in- 
creases in the price of energy carriers unlock additional 
energy resources by enabling a higher marginal cost of 
production to become economical. While low cost stock 
resources become increasingly depleted, increases in price 
of energy carriers enable the exploitation of additional high 



cost resources in order to supply the demand. It may thus 
be inferred that price paths in time produce supply paths 
in time for energy resources. It is however the demand for 
energy carriers that determine their prices: they adjust 
in such a way that the supply resulting from the sum of 
resources unlocked meets the price adjusted demand. 8 

Following the terminology of Mckelvey (1972) and 
Rogncr (1997), reserves are seen as continuously con- 
sumed but also continuously expanded at the expense of 
resources. Reserves are defined as those currently eco- 
nomical to exploit, and the boundary delimiting reserves 
from resources is defined by prices, which evolves in order 
to maintain a certain level of reserves out of existing re- 
sources, with respect to the demand. As we show below, 
it is remarkable that on the global level, the ratio between 
the rate of consumption of oil and gas to the size of their 
associated reserves has been nearly constant in recent his- 
tory. This aspect strongly strengthens the assertion that it 
is the evolution of prices that enables the size of reserves 
to follow the magnitude of their respective consumption 
levels, which we take as a starting postulate in order to 
define the theory given here. 

Meanwhile, reserves can expand with discoveries and 
technological change. Discoveries of resources that are not 
in areas where they are expected to be found (i.e. not 
counted as inferred or speculative resources) are part of the 
unknown/unknowable resource endowment, and therefore 
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may be considered as uncertainty. Similarly, technological 
breakthroughs in the oil industry are not straightforward 
to predict, and are also considered part of the uncertainty. 
A complete description of our approach to uncertainty is 
given in Mercure and Salas (2012). A Monte-Carlo simu- 
lations approach may be used by defining large numbers 
of variations of cost distributions of resources actually in- 
place, whether discovered or not, and generating differ- 
ent scenarios of resource use for different levels of scarcity. 
Therefore, in the model presented here, the process of dis- 
covery is not assumed to take an important role, or, in 
other words, the process of discovery is considered part of 
the gradual resource consumption process. 

2.2. Calculation of energy flows from existing stocks 

Flows of stock resources and associated depletion can 
be calculated given time paths for their associated carrier 
prices which unlock just the right amount of energy at 
every time step that is not already produced by other re- 
sources types, in order to supply the total energy demand. 
This increase in price is associated with the marginal cost 
of production for this energy resource. 

Following the first panel of figure 1 (see also figure 1 in 
Mercure and Salas, 2012), we take uq(C) as the initial cost 
distribution of a particular type of resource (e.g. oil, gas, 
coal or uranium) , and a time dependent resource distribu- 
tion n(C, t) which represents the cost distributed amounts 
of resource left at time t. n(C,t) is equivalent to a his- 
togram of all units (e.g. barrels, tons, etc) of a particular 
resource type that are assumed to be in place and ranked 
according to their cost of exploitation. We take the as- 
sumption that the rate of extraction of stock resources in 
each cost range (i.e. between C and C + dC) , at any time 
and price, is proportional to the amount of resource left in 
that cost range, if the latter is considered economical to 
extract, with a probability that it is considered so. 9 We 
take a continuous step-like function f(P(t) — C), 10 which 
equals one below this boundary and zero above, as the 
probability of extraction of resources in the cost between 
C and C + dC . This is shown in panel a) of figure 1, 
where a hypothetical resource distribution is shown as a 
solid black line, f(P(t) — C) is shown as a red curve that 
converges to one towards low values of C, and reserves cor- 
respond to the section of the distribution situated to the 
left of this curve. 

While the remaining cost distributed resources 
are n(C, t), 11 cost distributed reserves correspond to 
n(C, t)f(P(t) — C). If the constant fraction vq of reserves 
in each cost range are extracted at any time, the flow of 
resources is therefore as follows: 



dn(C,t) 
dt 



-v n(C,t)f(P(t)-C), 



The time dependent energy carrier price P(t) drives the 
evolution in cost of the boundary between reserves and re- 
sources. For a constant or increasing commodity demand, 



as the size of reserves decreases following consumption, 
this flow decreases and generates an upwards movement 
of P(t), shown in panel b), where the distribution of re- 
sources decreases in the low cost range and the boundary 
f(P(t) — C) moves to the right. This produces a time 
dependent supply (or flow) F{t) which is the sum of the 
production in all extraction cost ranges during the unit 
time dt: 



F(t) = 



dN(t) 
df~ 



dn(C,t) 
dt 



da 



(2) 



providing a connection between commodity prices P(t) 
and commodity flows F{t). Note however that this connec- 
tion does not have a unique functional form, but changes 
depending on the history of P(t) and the amount of re- 
sources left, a fact due to the integral of eq. 2, where for 
instance completely different historical paths and values of 
P at a particular time can lead to the same value of F. 
This indicates possible hysteresis and corresponds to path 
dependence. Details of these properties arc given in the 
supplementary information. 

These equations can be used in two ways, depending 
which variable is taken as exogenous and which is endoge- 
nous. For an assumed time dependent price P(t), the flow 
F(t) is straightforward to calculate numerically using a 
discrete time step using equations 1 and 2 (the forward 
problem). This is shown in panel c), where a range of 
time dependent flows is depicted, associated with the un- 
certainty in the amount of resources that actually turn 
out to be in place, using a linearly increasing commodity 
price. This flow is low at low price values, where no re- 
sources are economical to exploit, and low at high prices, 
where all resources have been consumed. 

Conversely, for an assumed commodity demand, the 
price value that unlocks just the right amount of resources 
can be obtained by trial and error using an optimisa- 
tion procedure (the reverse problem). 12 Panel d) depicts 
schematically the result of such a calculation with a range 
associated with the uncertain amount of resources in place, 
where lower prices result from higher amounts resources 
and the reverse, using as an example a large constant 
commodity demand. This results in a price value that 
gradually increases as resources are consumed, accelerat- 
ing when remaining resources are small compared to the 
level of consumption, forcing the exploitation of expensive 
resource occurrences. It eventually diverges as depletion 
sets in. In both examples, the red curves delimit the 96% 
confidence region associated with resource assessment un- 
certainty, while the blue curve denotes the most probable 
values. 



(1) 2.3. Constant global production to reserve ratios 



Insight for justifying the theory presented above can be 
obtained by considering the size of reserves for energy re- 
sources and their rates of exploitation. Reserves continu- 
ously expand, but are continuously consumed as well (for 
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Figure 2: Reserve to production ratios for stock resources for Oil (Left) and Gas (Right). This ratio corresponds to the inverse of uo, in 
years. In the case of oil, this converges globally towards 42 years after around 1987, and is constant at 62 years for gas throughout the period. 



a discussion see Mckelvey, 1972). The size of global oil 
and gas reserves has never been constant in history, nei- 
ther has global production. However, the ratios of pro- 
duction to reserves, when taken at the global level (only), 
are remarkably constant, given in figure 2, for which his- 
torical reserve and production data for oil and gas were 
obtained from the BP Statistical Review of World Energy 
(BP, 2009). 

Reserves concern a particular subset of specific resources 
(oil, gas, etc) that is both very probable to be in place and 
feature costs of extraction that are within a range deemed 
competitive given the price of internationally traded com- 
modities. The ratio of the rate at which reserves are con- 
sumed to the size of the reserves is equal to vq, the pa- 
rameter that controls the behaviour of equation 1 (details 
on the behaviour of the equation for different values of 
Uq are given in the supplementary information). On the 
regional level, various non-constant values for i/q are ob- 
served (sec fig. 2). Different regions have different energy 
policies related to their own political and geophysical sit- 
uations. However, their aggregated output supplies the 
international demand for resources, which sets the price. 
This price moves up and down in order to adjust the upper 
value of cost that enables resources to be extracted; it de- 
fines the size of reserves out of the resource base. A large 
amount of trade occurs between regions of the world, and 
overall the value of v$, on the global level, remains remark- 
ably constant. This supports the fact that the perceived 
price level, excluding short term fluctuations from specu- 
lation and hoarding except for a cost associated with risk 
due to volatility, evolves such that the ratio between pro- 
duction and reserves remains approximately constant, by 
unlocking just the right amount of resources to include into 
reserves in order to supply demand, given that reserves are 
exploited at a rate of z^o- 13 It also reflects empirically the 
maximum rate at which reserves can be extracted. 

In the case of oil, (left panel of fig. 2), Vq converges 



towards 42 years and becomes constant after 1987. Its 
evolution before this date is related to the oil shocks of the 
late 1970s, where OPEC was formed and Middle Eastern 
production decreased faster than reserves. In the case of 
gas, (right panel of fig. 2), this ratio is constant at 62 years 
throughout the period. These values correspond to the 
time it would take to deplete reserves if the price did not 
increase in order to unlock additional resources. In the 
case of coal and uranium, no historical reserve data is given 
by BP, but a similar process is likely to take place, and 
for which the 2008 values are of Vq 1 of 122 and 16 years 
respectively. 14 These values are assumed to remain true 
for the future in the calculations that follow. 

The significance of this phenomenon is not as surprising 
as it may initially appear. BP reserve data corresponds 
to a perception by the industry of the current outlook of 
energy resources, and their expectations regarding prices 
and global demand. It is natural to expect energy firms to 
increasingly expand their own conception of the economi- 
cal limit of reserves by considering to develop new projects 
that were considered prohibitively expensive in the past, 
in order to maintain reserves to a certain level, and that 
this level should evolve following global demand. This in- 
cludes for instance the current arctic oil exploration frenzy, 
Canadian tar sands, ultra-deep offshore rigs, etc. 

3. Global overview of stock energy resources 

3.1. Assumptions 

In this section, two extreme modelling exercises are car- 
ried out to explore the model properties before provid- 
ing an example of projection produced within a broader 
modelling framework. The assumptions for these exercises 
and for the projections are given here. The exploration of 
possible of stock resource flows and associated world mar- 
kets can be performed using either prices as exogenous in 
order to calculate resulting flows, or using flows (energy 
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Figure 3: Assumptions for calculations of stock resource flows and prices. Left Exogenous linear marginal cost component of price assumptions 
used to calculate flows, right Exogenous energy demand components used to calculate marginal cost components of prices. Note that here, 
the nuclear contribution to the total primary energy demand (TPED) appears as 2% instead as the 6% quoted in IEA TPED data IEA 
(2010) (IEA data for nuclear electricity production is 10 EJ). This is due to IEA's use of an arbitrary efficiency factor for converting 10 EJ 
of nuclear electricity into 30 EJ of primary nuclear fuel, out of 514 EJ of TPED. The value used here is 10 EJ out of 494 EJ of TPED, since 
the conversion factor of nuclear reactors is already included in the resource data of Part I. 



demand) as exogenous in order to calculate commodity 
prices, excluding in both cases the effects of hoarding and 
short term demand fluctuations. Values used or produced 
here and henceforth in this work correspond to marginal 
costs of production, rather than real prices, and whatever 
margins of profit, fixed costs and other cost component as 
well as fluctuations may be added to these marginal costs 
in order to construct real endogenous prices (prices cal- 
culated within an energy model as opposed to assumed), 
using separate assumptions, not done here, but left to the 
discretion of the modeller. 

As a first exercise, we present in section 3.2 a calculation 
of stock resource flows given assumptions for the marginal 
costs deemed economic given energy carrier prices, denoted 
P. The assumptions are given in the left panel of fig- 
ure 3. Values for vq are given in the preceding section. 
Given these, starting marginal costs required for supply- 
ing current demand were evaluated by finding which value 
of P generates a dN/dt that equals the current demand. 
Following this, rates for the increase of the carrier prices 
were used that generate an increase in supply consistent 
with current total primary energy demand increase, but 
were maintained constant throughout the century (linear 
prices). The resulting flows are given. 

As a second exercise, we present in section 3.3 a calcula- 
tion of the required marginal costs of energy production for 
each stock resource for a scenario where total primary en- 
ergy demand increases to around 900 EJ/y, but where the 
current shares of energy demand for these resources (coal, 
oil, gas, uranium) are rigidly maintained until 2100 (i.e. 
the structure of the current energy system is maintained) . 
This energy demand curve until 2100 is within a range con- 
sistent with many recent projections, a review of which is 



given by Edenhofcr et al. (2010). The results are given in 
section 3.3. These rigid assumptions specifically exclude 
technology substitutions in both sets of calculations, an 
aspect which is explored separately in section 3.4. 

Neither of these calculations produce realistic scenar- 
ios of energy use and prices, as we demonstrate in sec- 
tion 3.4. This is due to the absence of technology sub- 
stitutions that enable to avoid price escalations related to 
scarcity by switching away from these sources, and of a 
dynamic feedback with energy demand from the economy. 
While the second aspect cannot be analysed here with- 
out a full blown description of extensive global macroe- 
conomic modelling (e.g. using E3MG), the first is readily 
explored, however in the power sector only, by introducing 
the model described above into FTT:Power. FTT:Power 
enables to explore technology substitutions that are likely 
to arise in the power sector following changes in relative 
costs of energy systems, which include the cost of non- 
renewable fuels. Therefore, in the event of depletion and 
scarcity of some resources, the model endogenously gener- 
ates switches of technology by gradually phasing out some 
types of systems, following possible rates of diffusion of 
new technology and rates of decommission of old systems 
(see Mercure, 2012a, for a complete description of the 
model). Such changes reduce the demand for some types 
of resources, avoiding their depletion and large price in- 
creases that would occur in the case of a perfectly rigid 
demand. The interaction of FTT: Power with E3MG will 
generate complex energy-economy interactions which will 
be explored in detail elsewhere. Thus, for this work and 
for simplicity, two global energy demand scenarios were 
chosen based on IEA projections. The first is a baseline 
scenario of global energy policy grounded on the assump- 
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tions of the IEA's World Energy Outlook new policies sce- 
nario (see IEA, 2010, appendix B, for a detailed descrip- 
tion). The second corresponds to a strong mitigation sce- 
nario aiming at reducing GHG emissions to below 50% 
of the 1990 levels, following the assumptions of the 450 
scenario of the World Energy Outlook, however with ad- 
ditional regulations and support for particular renewable 
technologies. Out of these scenarios, all exogenous inputs 
to FTT:Power were taken. This corresponds to the global 
demand for electricity and the global demand for energy 
resources that do not originate from the power sector (i.e. 
final use of coal, oil and gas by industries, households and 
transport). These assumptions are given in section 3.4. 

3.2. Flows from stock energy sources for exogenous prices 

In this first extreme modelling exercise, scenarios for 
flows of stock resources were produced for oil, gas, coal 
and uranium using the resource distributions underlying 
the cost-supply curves given in figure 4 of our previous 
work (Mercure and Salas, 2012), current rates of resource 
exploitation (values for vq given above) and linear extrap- 
olations of energy carrier marginal costs of production. In 
the case of coal, the sum of the distributions for hard and 
soft coal was used. The results are given in figure 4. In 
all cases, the curves start in 2008 at the current global 
production values reported by the IEA (2010), 15 with zero 
uncertainty (i.e. current reserves are known). As time 
progresses, the increasing uncertainty in resources assess- 
ments at higher levels of use produces ever larger ranges of 
possible resource production values, or ranges of possible 
consumption paths, delimited by the red curves. These 
however must ultimately converge back to low flow values 
when peaking and depletion occurs. 

Different results are obtained depending on the size of 
the various stocks. In the case of oil, which includes all 
types of unconventional oil, a peak in production occurs 
at around 2060, after which depletion begins. A similar 
situation occurs with natural gas, which includes all types 
of unconventional gas and methane hydrates, peaking later 
near 2080. Coal resources, however, are very large and 
depletion does not occur within a foreseeable future. It 
can only do so very far outside the time horizon of 2100. 
Resources of natural uranium, as reported by the IAEA 
(2009), are found to become depleted rapidly within the 
current century after peaking before 2025. 16 

Potential flow values vary highly between resource types. 
Projected flows from oil resources are the largest, up to 
600 EJ/y, giving however a faster rate of depletion com- 
pared to coal and natural gas. This is due to the current 
large rate of extraction to resource ratio vq . Their massive 
expansion occurring after 2020 is due to the large scale ex- 
ploitation of unconventional oil such as the tar sands and 
oil shales. While natural gas resources are smaller than 
those of oil, their depletion is projected farther into the 
future due to a lower extraction rate to resource ratio vq. 



Their massive expansion after 2040 is related to large scale 
exploitation of unconventional sources such as shale gas. 

In the case of coal, the rate of extraction is similar to 
that of gas, but their reserves are much larger, projecting 
the depletion far beyond 2100. In the case of uranium, 
given the small resource base, the low burn-up rates of 
current thermal reactors and the high value for uq, the ex- 
pected flow is very small compared to those of fossil fuels, 
but the depletion is projected to occur rapidly within this 
century This indicates that dramatically higher conver- 
sion efficiencies are necessary to extend the resource base 
beyond the end of this century, which can be achieved with 
fast breeder reactors (Bonche, 2002, Nuttall, 2005). 

Note that these flow paths are not forecasts in any way. 
They are possible scenarios of resource use, given known 
(and uncertain) resources bases and realisable extraction 
rates and price paths, even though the price of carriers is 
hardly likely to follow a linear trend. Nevertheless, for any 
resource extraction path, the cumulative flow up to com- 
plete depletion must be equal to the technical potential of 
the resource, a requirement that has been carefully verified 
in these calculations. Therefore, for higher rates of extrac- 
tion and resulting higher resource flows, depiction must 
occur slightly sooner than depicted here, and conversely, 
for slower resource use, depletion may occur slightly later. 

3. 3. Price paths for exogenous flows 

Conversely to the previous section, as a second extreme 
modelling exercise, the reverse problem is posed where one 
looks for the appropriate price of resources that unlocks 
just the right amount of resources to meet an exogenous 
(rigid 17 ) demand. In this case, the demand, or resource 
flow, is given as exogenous and the price, or marginal cost, 
is evaluated. This is done by performing an inverse calcula- 
tion using a price optimisation of eq. 1, such that the value 
of eq. 2 is equal to the demand, separately for each type 
of energy carrier (oil, gas, coal, U). In such scenarios, it is 
possible that at a certain point in time, given the values 
of fo, the remaining resource base cannot meet anymore 
the demand. In such a case, the price values gradually run 
away to infinity, signifying that all resources situated at 
all possible prices of extraction are under intense exploita- 
tion. Such a situation is very unlikely to occur, since the 
opportunity cost of using these expensive resources would 
become very large and other technologies and energy re- 
sources would be more cost effective, leading to technology 
substitution, switching away from that resource before it 
runs out. Alternatively, the global economy may also read- 
just its energy demand in order to avoid diverging prices 
of energy commodities. The divergence of prices therefore 
stems from rigid commodity demand values that do not 
respond to price signals. 

Figure 5 presents the results of such an exercise, us- 
ing the assumptions described in the right panel of fig. 3, 
where the current composition of the energy sector is main- 
tained with a total energy demand scaling up to near 
900 EJ/y in 2100. The blue curves correspond to the most 
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Figure 4: Flows of oil, natural gas, coal and uranium, as calculated using equation 1 and the cost-supply curves given in figure 4 in 
Mercure and Salas (2012). The 95% confidence level region is situated between the red curves, while the blue curves indicate the most 
probable flow values. The curves start in 2008 at current energy consumption values, given by the IEA. 



probable resource bases given by the blue curves in fig. 4 
of Mercure and Salas (2012). Meanwhile, the red curves 
delimit the 96% confidence region, where the upper red 
curves correspond to the lower bounds for resources, while 
the lower red curves correspond to the higher end of the 
resource ranges. Therefore, in all cases, the marginal cost 
values calculated in the low end of the resource ranges di- 
verge, while most of the curves for the upper ranges do 
not. 

In the case of oil and gas, gradual increases are observed 
in the marginal cost values, with a change in slope occur- 
ring between 2020 and 2030. This is related to the price 
enabling the accession of large unconventional resources, 
which include predominantly oil sands and shale gas re- 
spectively. The availability of these large resources tend 
to damp out possible future increases in price. 18 Mean- 



while, in the case of coal, the marginal cost value is hardly 
affected by demand at all, unless resources turn out much 
smaller than expected, which is very unlikely. This is due 
to the very large resource being situated in a narrow range 
of extraction costs. 19 

Finally, in the case of uranium, the marginal cost is ex- 
pected to diverge, and the resource to run out, over the 
whole resource uncertainty range, if the current share of 
energy demand supplied by nuclear is maintained up to 
2100, using current technology without recycling waste. At 
the current uranium burn-up rates, the resource base is in- 
sufficient. This indicates that the nuclear industry should 
cither decrease its share of electricity generation signifi- 
cantly by 2100, or that much higher efficiency rates in 
resource used per unit of electricity produced are achieved 
before then, involving possibly a much higher rate of re- 
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Figure 5: Marginal cost calculations performed for four stock energy resources using the assumptions described in the right panel of fig. 3, 
which maintains the current composition of the energy sector, while expanding the total energy demand up to nearly 1000 EJ, using the 
resource data from Mercure and Salas (2012). The blue curves correspond to the most probable values given resource uncertainty, while the 
red curves delimit the 96% confidence range, where lower amounts lead to marginal costs increasing earlier. 



cycling of nuclear waste than occurs at present. Future 
energy systems planners will adopt either of these solu- 
tions in order to avoid this projected fuel cost escalation. 

3.4- Real systems: allowing technology substitutions 

Although the modelling exercises given in the last two 
sections provide insight on the scale of available resources 
and on the process of their gradual consumption, they both 
depict limiting situations that are very unlikely to occur. 
This is due to the facts that: 

1- There exists a feedback between prices and demand in 

the global economy 

2- Technology/resource substitution processes occur that 

enable reductions of the demand for specific commodi- 
ties. 

Prices are not likely to remain strictly linearly increasing 
as in section 3.2, and the composition of the demand is not 



likely to remain fixed in the future, as in section 3.3. Effec- 
tively, as depletion progresses, marginal costs of exploita- 
tion increase and prices increase, and these induce gradual 
technology switching and/or reductions in overall energy 
demand. Complete technology switching away from a par- 
ticular fuel occurs when the price of this fuel makes its use 
uneconomical. Therefore, the prices can never escalate to 
very high values as long as technology switching options 
exist since switching away occurs before the price diverges 
at complete depletion, and thus stock resources are never 
depicted entirely. Technology switching is however con- 
strained by capital lifetimes and can therefore take some 
time to take place, in particular in the power sector. 

The inverse calculation problem described in section 2.2 
to derive a marginal cost of production given an exogenous 
demand was introduced into the model of the global power 
sector FTT:Powcr, which specifically simulates technology 
switching given plant lifetimes and dynamic rates of tech- 
nology diffusion, described in length in Mercure (2012a). 
It thus provides an appropriate testing ground for this the- 
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Figure 6: Example of changes in the power sector which could generate strong increases in the marginal cost of natural gas and uranium, in 
comparison to a baseline scenario, where in the top panels electricity generation is given, while in the bottom panels the associated emissions 
are shown. On the left panels is shown the baseline scenario, which consists primarily in assuming that current policies extend into the 
future. On the right is shown a mitigation scenario where carbon pricing exists (starting at 22$/tC02 and increasing by 2% per year ) and 
support for the wind industry and CCS technology. Dashed black curves correspond to the same baseline as in the left panels, given for 
comparison. The 2008 technology emission factors were obtained by comparing emissions to electricity generation in IEA data, maintained 
in the projection. Horizontal dashed lines give reference to the 1990 levels, while vertical dashed separate projections from historical data. 



ory, an exercise that also generates direct insight on the 
effect to prices of the future composition of the power sys- 
tem constrained by natural resources. Figures 6 and 7 
present the results of simulations performed with this ver- 
sion of FTT: Power for all four non-renewable resources, 
oil, coal, gas and uranium, and the values of evaluated 
above. The particular example given was chosen on the 
basis that it describes well the properties of this model in 
order to demonstrate its validity, not for the goal of recom- 
mending any particular power sector technology scenario. 

The theory presented here however requires the global 
demand from all sectors for natural resources, not entirely 
provided by FTT:Power, which treats the power sector 
only. Thus, the missing components of the demand unre- 
lated to electricity production had to be taken as assump- 
tions, described below. 20 In all other respects, the model 
assumptions are very similar to those of the baseline pre- 
sented in Mercure (2012a). 21 These simulations were per- 
formed for two sets of policy and demand assumptions, 
a baseline scenario and a mitigation scenario. The top 
panels of figure 6 present electricity generation by type 
of technology for both scenarios, the baseline on the left 
and the mitigation scenario on the right. In these plots, 
the dashed vertical lines provide a visual reference to the 



present and the data to the left of these lines corresponds 
to historical data, while on the right are given the calcu- 
lated projections. Meanwhile, in the bottom panels, the 
associated emissions from fuel combustion in the power 
sector are given, where the horizontal dashed line indicates 
the 1990 level. 

The demand for energy commodities (coal, oil, gas and 
uranium) was calculated using endogenous values for fuel 
demand by the power sector, and exogenous values for fuel 
demand from the rest of the economy. Demand values for 
the rest of the economy are given in the top panels of 
figure 7 for each scenario. In the baseline, the demand 
for oil, originating primarily in the transport sector, was 
assumed to peak late in the century, motivated by a grad- 
ual transition to alternative transport technologies. The 
non power-related demand for gas, originating primarily 
from the industrial and buildings sectors, was assumed 
to increase gradually up to 2100, although slowing down 
due to gradually increasing overall efficiency in parallel 
with an increasing demand for heating services. The non- 
power demand for coal, originating primarily in the indus- 
trial sector, was assumed to rapidly peak and gradually 
decrease due to technology switching and increased use 
of natural gas. Finally, the demand for biomass was as- 
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Figure 7: Top panels exogenous assumptions for the demand for coal, gas, oil, uranium and biomass, that does not originate from the power 
sector, for the baseline (left) and the mitigation (right) scenarios given in figure 6. Bottom panels Marginal cost of production for coal, gas, 
oil and uranium that result from the global demand and the theory presented above, along with the marginal cost of biomass production 
which is obtained directly from the cost-supply curve. 



sumed to gradually increase, at a rate accelerating in the 
second half of the century due to a higher diffusion of bio- 
fuels for transport. In the mitigation scenario, oil demand 
peaks rapidly due to massive technology switching in the 
transport sector towards biofuel and electric cars. The 
demand for biomass increases sharply to supply this ad- 
ditional demand. The non-power demand for natural gas 
however remains relatively constant and the demand for 
coal sharply declines due to massive technology switching 
and electrification of industrial processes. 

The bottom panels of figure 7 present the resulting 
marginal cost of production of coal, gas, oil and uranium 
in these model runs of FTT:Power including the present 
theory for both scenarios, which can be used, with addi- 
tional chosen assumptions, to construct endogenous prices. 
In the baseline scenario, a strong increase in the cost of 
natural gas is observed, associated to an increasing global 
demand, forcing the price to enable the extraction of shale 
gas and more expensive resources. This rate of increase is 
however damped past 2020-2030 due to the large amount 
of resource available at these cost ranges. Meanwhile, the 
marginal cost of coal hardly changes, irrespective of the 
sharply increasing demand, reflecting the sheer scale of low 
cost coal resources. Nuclear reactors are mostly decom- 
missioned and see a sharp decline in the baseline scenario, 
leading to a decrease in the price of uranium. Finally, 



the price of oil increases due to increasing depletion, but 
the rate of increase is damped by the accession to mas- 
sive amounts of unconventional oil. However, since the 
main component of the demand does not originate from 
the power sector but the exogenous transport demand, the 
analysis of oil demand is outside the scope of this work. 22 

In the mitigation scenario, strong support is given to 
wind energy through a subsidy (35% of the LCOE through- 
out the simulation period), as well as through the pricing 
of CO2 emissions (starting at 22 $/t and increasing by 2% 
per year up to a value of 140 $/t in 2100), while mod- 
erate support is given to electricity production using cap- 
ture and storage (CCS) technology (10% of their respective 
LCOEs). 23 The introduction of large amounts of variable 
renewable electricity into the grid requires increases in the 
amount of flexible type of generation, which can be pro- 
vided for by, for instance, gas turbines, oil plants or hydro- 
electricity This motivates a massive expansion of the gas 
turbine technology into the electricity market, which even- 
tually dominates. The carbon pricing however motivates 
the installation of CCS on all gas turbines by 2060, reduc- 
ing drastically emissions. Meanwhile, the pricing of car- 
bon makes coal technologies gradually come out of favour, 
while the nuclear industry maintains a constant market 
share. However, with the assumption that nuclear reac- 
tors only use natural uranium, do not recycle waste and 
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maintain the low conversion efficiency of ordinary thermal 
reactors, the resource base is seen unable to maintain the 
share of nuclear capacity and a strong price increase for 
natural uranium is observed, which generates a decline of 
the nuclear industry starting at around 20 70. 24 The mas- 
sive expansion of gas turbines generates a stronger increase 
in the cost of gas compared to the baseline, which is how- 
ever damped slightly due to the large amounts of shale gas 
available. The cost of coal resources gradually decreases 
following the decline of coal in electricity generation. With 
a smaller demand for oil resources by the transport sector 
as assumed exogenously, the price of oil initially increases 
but stabilises and decreases slightly in the middle of the 
century. 

Finally, global power sector emissions are given for both 
scenarios in the bottom panels of figure 6. While emis- 
sions increase monotonically up to 37 Gt/y in 2100 in the 
baseline scenario, they peak at 14 Gt/y in the mitigation 
scenario in around 2030, they decrease afterwards mono- 
tonically and reach the 1990 level of 7 Gt in 2060, and then 
decrease to a low level of 1 Gt in 2080 where it remains con- 
stant. Thus the policy assumptions for the mitigation sce- 
nario are not stringent enough to reach 50% reductions by 
2050. Note that the technology mix outcome of the various 
possible subsidy schemes transform this situation greatly. 
While this example is not particularly attractive for pol- 
icy, it was chosen because it displays strong marginal cost 
changes related to changes in demand, where for instance, 
the development of biomass based electricity was pushed 
out of the market by the strong mutually beneficial com- 
bination of wind and gas turbines, which is not necessarily 
beneficial for overall emissions. 25 Detailed scenario analy- 
ses using FTT:Power are not the primary objective of this 
work and will be presented elsewhere. This scenario does, 
however, warn of the potentially strong effect that energy 
policy and future demand can have on the price of energy 
commodities. 

This section demonstrates that the combination of the 
theory presented above for treating dynamically the con- 
sumption of stock resources, with a model of technology 
substitution, enables to effectively project future marginal 
costs of energy commodities given exogenous demand val- 
ues from the rest of the economy. However, since technol- 
ogy substitutions can occur in all sectors of the economy, 
additional flexibility in demand values exists in the global 
economy. Therefore, while the power sector amounts to a 
very significant share of global fuel use, more accurate cal- 
culations for energy commodity prices can be performed 
using a combination of this theory with a complete fam- 
ily of models of technology substitution for all major fuel 
users of the global economy, for instance the FTT family 
(for the theoretical framework underlying the FTT family, 
see Mercure, 2012b), the development of which will enable 
to remove one by one the exogenous demand assumptions 
given above. This is a substantial project which is cur- 
rently under way and will be the subject of forthcoming 
publications. 



3.5. Impacts on climate policy 

Carefully designed climate policy must take into con- 
sideration the amounts of low emission energy resources 
available in the world in order to produce its desired out- 
comes. Additionally to this, however, climate policy mak- 
ing must evaluate its own effects on energy prices, the ef- 
fects of energy prices on technology substitution and onto 
the economy. 

Four aspects must be considered in such policy frame- 
works: 

1- The cost distributions of energy resources, 

2- The relationship between energy prices, the size of the 

resource bases and rates of exploitation, 

3- The damping effect of technology substitution on these 

prices, 

4- The effects of climate policy onto energy prices and 
their subsequent consequences onto the well-being of 
the global economy, including through the price of 
electricity. 

The theory presented above, in combination with our pre- 
vious work (Mercure, 2012a,b, Mercure and Salas, 2012) 
is appropriate to treat these issues, for instance by using 
the model FTT:Powcr, excluding the effects of high en- 
ergy prices on the economy, which can be modelled with 
an economic model such as E3MG. 

Failing to address these issues, by for instance being op- 
timistic on the amounts of low carbon resources available, 
is likely to lead to badly planned energy policy where high 
energy carrier prices result (e.g. the price of uranium) 
and/or strong rebound effects arise due to reductions in 
prices, generating additional demand for these carriers if 
the price difference is not absorbed by energy pricing pol- 
icy (e.g. energy taxes), where for instance reducing the 
consumption of oil in one region of the world reduces the oil 
price, and induces additional consumption elsewhere. Op- 
timistic reports on the availability of low carbon resources 
abound in the literature, including the World Energy As- 
sessment (UNDP, 2000) and the IPCC special report onre- 
newables (IPCC, 2011), which list energy economic poten- 
tials as single values rather than cost distributions, making 
these considerations ambiguous. Additionally, reports ex- 
ist that claim the feasibility of the complete replacement 
of fossil fuels by renewables by 2050 including a report 
by the World Wildlife Fund (WWF, 2011), which do not 
take consideration of either cost distributions of resources, 
provide reliable estimates of energy resources or treat the 
feasibility of massive rapid diffusion of new technology into 
the marketplace. Such overly optimistic studies must be 
taken with care by the policy community. Energy systems 
are complex and deserve to be treated using complex dy- 
namic modelling. 
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4. Conclusion 

We have presented a model that, along with our previous 
paper (Mercure and Salas, 2012), completes our work on 
the economic potentials of energy resources by generating 
marginal cost of non-renewable energy carrier production, 
which can be used to project future energy prices. We 
have defined a model, based on cost distributions of non- 
renewable energy resources published in previous work, 
that simulates the simultaneous consumption and expan- 
sion of energy reserves. This model uses a dynamic differ- 
ential equation approach to calculate at each time step of 
a projection the amounts of energy resources consumed at 
a particular price. This generates an endogenous resource 
consumption path given an exogenous price path. This 
model can however be used in reverse, where the price that 
generates an appropriate resource supply is found by opti- 
misation, generating an endogenous price calculation given 
an exogenous demand. Using both approaches, ranges of 
possible consumption paths as well as possible marginal 
cost ranges have been derived given the available resource 
bases and associated uncertainty ranges. The use of rigid 
demand or price assumptions however lead to unrealis- 
tic price escalations or depletion timescales. By connect- 
ing this model to one of technology substitution based on 
rules of investor choice over pairwise comparisons, price 
increases are damped through substitution, avoiding the 
depletion of particular resources and the associated price 
increases. The analysis presented thus generates insight 
for energy planning related to non-renewable resources. It 
generates consumption peaking timescales for oil, natural 
gas, coal and uranium resources, where in particular, em- 
phasis is given to the limited amounts of available natural 
uranium. When enabling technological change, technol- 
ogy substitutions are observed that avoid marginal cost 
increases. Reflections on the policy implications of the 
model and results presented are provided. 

This model forms the third part of a theoretical frame- 
work for exploring global future technology transforma- 
tions and climate policy, which consists in (1) a technol- 
ogy substitution model based on a coupled family of logis- 
tic differential equations representing technology diffusion 
into the marketplace, (2) a database and theoretical ba- 
sis for tracking the use of global energy resources, and 
(3) this work, which provides the basis for a model for 
endogenous price formation. This cornerstone will enable 
the completion of the full FTT family of technology models 
which aims at projecting multi-sectoral global greenhouse 
gas emissions. 
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Notes 

1 The global output of the energy and fuel supply industries makes 
2-3% of global GDP and decreasing, values obtained from E3MG 
calculations by T. Barker (2012), private communication. 

2 This is also a pronounced effect in E3MG results. 

3 Most models rely on outdated and fixed (i.e. not time depen- 
dent) cost-supply curves from Rogner (1997). 

4 As can readily be observed using E3MG with different scenarios 
of energy prices. E3MG is an econometric model that extrapolates 
such trends from data. 

5 Power sector technology substitution can occur at timescales not 
much shorter than 30-40 years. 

6 This can be seen with a dip in oil sands production at the dip in 
crude oil price in IEA oil price and production data around 2008-2009 
(IEA, 2011). 

7 It does not replicate the reasons why people hoard stockpiles or 
the expectations of political conflicts or supply restrictions related to 

■/^iWA0^mRd^i&iim^ a ISkf}^d9 vcr a medium term artificial 
demand component (positive or negative) originating from changes 
in stockpiles, the model can be used to predict medium term price 
fluctuations, but not price hikes associated to very short term supply 
problems. 

8 Strictly speaking, this is true over the medium term where de- 
mand adjustments have time to take place. 

9 For example, the supply is proportional to the number and/or 
size of wells or mines, which is proportional to the amount of eco- 
nomical resources (reserves) left. 

10 The function can be a smooth rounded step to reflect variations 
over the response to the price, see the supplementary information for 
details. This is similar to particular conceptual problem in reaction 
chemistry and physics, see Mercure et al. (2005). 



16 



llr rhe total remaining resources correspond to the integral of n, 
f™n(C,t)dC. 

12 This is done by trial and error because it is not possible to know 
the integrand of eq. 2, which is function of C, from its integrated 
result F(t). One must therefore take guesses over which value of P(t), 
f(P(t) — C) and therefore of the integrand, that gives a particular 
value for F{t). 

13 vq can be seen as the time it would take to consume reserves 
at the current consumption level, and therefore has units of inverse 
years. 

14 Data for uranium are taken from both BP (2009) and IAEA 
(2009). 

15 Note that the IEA's value for nuclear electricity production of 
10 EJ was used, not its reported value of primary nuclear fuel of 
30 EJ. Efficiency factors for thermal reactors have already been taken 
into consideration in the cost-supply curve for uranium. 

16 this excludes the reuse of fissile material available in nuclear 
waste, which would probably become economical before the com- 
plete depletion of natural uranium resources. After 2025, the 
amounts of available fissile material in the waste produced by pre- 
vious use of natural U will be large, and the high costs of recycling 
nuclear waste will eventually be equalled by the increasing costs of 
mining U ore which will become ever more difficult to reach. 

17 Rigid in the sense that the demand does not respond to prices. 

18 In a scenario where no additional environmental regulations pre- 
vent their exploitation, an obviously disputable assumption. In the 
event where such regulations arise (limiting fracking for instance, or 
regulations being instated in Canada regarding river and land con- 
tamination from tar sands processing), the situation may become 
similar to the lower resource ranges given by the lower (high prices) 
red curves. 

19 Note that for coal, it is not necessary to include low grade re- 
sources in high cost ranges, since normal grade resources are very 
large and unlikely to be consumed entirely within this century. 

20 The future construction of models of technology substitutions in 
the transport, industrial and household sectors will in time replace 
these fixed assumptions. 

21 The baseline here features no carbon pricing and no subsidies or 
taxes on any technology. 

22 A similar treatment of technology substitution in the transport 
sector would be required in order to produce a dynamic demand for 
oil, electricity and biofucls that responds to prices, which will be the 
subject of future work and a new model, FTT: Transport. 

23 Additional subsidies are given to biomass based electricity and 
solar technologies, of 35% and 50% respectively, without much im- 
pact in this particular scenario. These are effectively pushed out of 
the market by wind and gas turbines. 

24 Note that although the price of U goes off the scale of the lower 
right panel of figure 7, it does not diverge to infinity. Its scale of 
increase is related to the rate at which the nuclear industry can 
decommission its power stations in a scenario where waste recycling 
is not allowed. 

25 For instance, negative emissions from biomass gasification with 
carbon capture and storage can reduce emissions by much larger 
amounts. 
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Supplementary Material 



S.l. Mathematical properties of the resource flow equation 

5.1.1. Introduction 

We present in this supplementary material a mathematical digression that explores the mathematical properties of 
equations 1 and 2 of the main text, which establish a relationship between the price P(t) of an energy carrier derived 
from a particular type of non- renewable energy resource (e.g. coal, oil, gas, U) and its consumption, or flow, F(t). As 
stated in the main text, the relationship is not functional, i.e. F(t) cannot be written as a single valued function of 
P(t) or the reverse. The relationship is path dependent, and therefore depends on the history of the system and on its 
starting point. It thus features hysteresis. The model is slightly inspired from a physical model of heat generation in the 
relaxation of an amorphous network, as given by Mercure et al. (2005). This material also presents the different limiting 
behaviours of this system for various values of its parameter Uq. 

5.1. 2. Cost-supply curves 

As was introduced in Mercure (2012a) and used in Mercure and Salas (2012), we define an initial (present day) 
distribution of non- renewable (stock) resource 77.0(C) function of cost C. This is a histogram of the number of resource 
units between cost values of C and C + dC. This is a density function; the amount of resources available at costs between 
the values of C\ and C2 is 



An example of such a starting resource distribution is given in the left panel of figure 8. At the present day, given the 
cost of extraction of every resource unit, the total amount TV of resources available below a cost value of C is 



Associated to this is the inverse relationship, the marginal cost of the resource C{N) given that N units have already 
been exploited; this is the cost-quantity curve. This is shown in the right panel of figure 8. 

The amount N of stock resources available below cost C cannot be exploited instantaneously however, and furthermore, 
their owner might not be willing to extract and sell them at the current price of the associated commodity. Therefore, 
this amount N cannot be sold for a the price associated with the marginal cost of production C. If less resources are 
extracted at prices below C than the total amount available in this cost range, resources situated higher up along the 
distribution must be used. Thus the cost-supply curve framework is not appropriate to use for stock resources with an 
international market. A similar statement could be invoked for renewable resources, however it has much less impact 
and the cost-supply curve framework is much more appropriate there. This is discussed in section S.l. 8. 

S.l. 3. The probability of extraction 

Each resource unit has at every instant a probability of being extracted, and this probability depends on its probable 
individual cost of extraction C and the market price of the commodity P(t), /(C, P(t)). Its cost of extraction is 
uncertain; therefore a probability distribution exists for the value of its cost of extraction h(C). Additionally, the market 
price is stochastic and has a certain volatility, with a certain standard deviation and probability distribution g (P(t)). 
The probability / (C, P(i)) is related to h(C) and g(P) through the integral of their convolution, providing the 'rounded 
step-like' probability function depicted in figure 2 of the main text. If h and g are normal distributions, then / is an 
error function with width equal to the root of the sum of the squares of the widths of h and g. These assumptions enable 
to avoid the use of a sharp step function, which would introduce unwanted kinks into calculations. 

The probability of a resource unit of being extracted can moreover be expressed in terms of the difference between its 
most probable marginal cost of production and the mean value of the price, / (P(t) — C). 

S.l. 4- A differential equation for resource flows 

While the initial (present day) distribution of resources is denoted with no(C), the distribution of future amounts of 
resources left as they are gradually consumed is denoted as the time dependent function n(C,t). This function depends 
on time in two ways, by itself and through the value of the price P(i), and thus strictly speaking, should be written 
as n(C,t, P(t)). This property is the one that leads to path dependence, since the time derivative involves two terms, 
shown below. 




(3) 




(4) 
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Figure 8: Left Example of cost distribution for a specific natural resource nrj(C) as a function of exploitation cost C. Right Cost-quantity 
curve for that resource type. 

While the cost distributed amounts of resources left is n(C,t,P(t)), not all resources are exploited but only those 
which have a probability of being exploited, given by 

n(C,t,P(t))f(P(t)-C). (5) 

This corresponds to cost distributed reserves, and the size of the reserves depend on the price, expanding when the price 
increases. If a constant fraction vq of reserves arc consumed with the time interval dt, then the flow of resources during 
that interval is 2 

dn (C, t, P(t)) = -v n (C, t) f (P(t) - C) dt, (6) 

which is the main equation of the model. This equation has no complete analytical solution but can be evaluated 
numerically, which is done for instance in FTT:Power using a discrete time step. Note that the negative sign stems from 
the fact that n corresponds to cost distributed resources that are left at time t, and that the change in resources left is 
negative. The flow of resources F(t) is the time derivative of the total amount of resources left at all cost values N(C,t), 
which itself is the cost integral of the resource distribution n(C, t). 

F (P(t), i) = = - jT ^§^dC = J™ u n (C, t) f (P(t) - C) dC. (7) 

It is important to emphasise that the notation F(P(t),t) signifies that F is not a single function of the price P(t), or 
conversely, since an infinite number of distributions, or integrands of eq. 7, can produce the same value for F. 

S.1.5. Mathematical properties when the price is constant 

A few important properties of equation 6 may be derived from the simple situation where the current price of oil P, 
is independent of time. Eq. 6 is solved simply: 

dn 

— = -v f(P - C)dt, (8) 
n 

^n(C,pt)=n (C)e~ tv ^ p - c \ (9) 

where no is the initial distribution of the resource previously defined, which we take, for simplicity for now, as a constant 
extending to very large values. We observe from this result that for a constant price, the amount of resource left at each 
cost value decreases exponentially in time, more and more slowly at higher and higher values of C. 



2 The actual fraction of cost distributed reserves consumed is actually uaf{P(t) — C), a slight subtlety, which could become important is 
f(P(t) — C) is very 'rounded'. 
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Figure 9: Left Distribution of oil resources n as a function of extraction cost C at different times. Right Total oil extracted as a function of 
time. 



The behaviour of n as a function of C is not quite as simple as this, but may be calculated numerically using a simple 
functional form for /, 3 and is shown in the left panel of fig. 9 for different times after exploitation began. The times 
noted as t\ to te increase exponentially by a factor of 2 between each. Starting from the initial distribution no, a certain 
range of low cost units are first extracted, until they are depleted at time £3. From then on, units at higher costs start 
to be exploited and the mid-point of n starts to move towards higher values of Co up to t±. shows n after a very 
long time. We observe that for t — > 00, all units are eventually used up and n = for all values of Cq. This is due to a 
non- vanishing value of the probability distribution function f(P — C). 4 

Furthermore, by calculating the area underneath uq — n(Co,t), we find the time dependence of the total number of 
units extracted. Fig. 9, right panel, shows this value against ]n(l + t). We observe that before the beginning of depletion, 
which occurs at around time £3, the extraction is fast and linear against ln(l + t). After time £3, it slows down but is 
again perfectly linear against ln(l + 1). We conclude that even though for t — > 00, the system uses all existing resources, 
the progression of the extraction becomes exponentially slower and slower. At any practical time, the shape of ro(Co, C, t) 
actually corresponds qualitatively to that at time £3 in fig. 9, where low cost units have been consumed, and progresses 
very slowly. 

Thus this demonstrates that consumption occurs even when the price does not increase. The supply, however, most 
likely does not meet the demand, since it monotonically decreases. 

S.1.6. Consumption for a changing price 

In the general case where the market price of units varies in time, many types of behaviour can be observed. The 
important mechanism at work, however, is that as low cost units are depleted, in order to keep a supply of units which 
docs not plummet, the price of the resource must increase in time in order to make the extraction of more expensive 
units economical. There is thus a direct relation between supply and price, but it cannot be expressed by a simple curve. 
We derive here the general solution to eq. 6. 

The general solution to eq. 6 is as follows: 

/ ^r=-v / f(P(t')-C)dt', (10) 

Jn n JO 

which yields 

n(C,P(t),t) =no(C)e- u % t/(p(t ' ) - c)dt '. (11) 

Two time dependences are specifically denoted, that inherent to the price P(t), and that associated to the exponentially 
decreasing exploitation that occurs at constant C. It is thus clear that the number of units n(C,P{t),t) at any time 



3 A logistic function was used in this particular case; however an error function may be as appropriate. 

4 The non-zero value of / concerns the non-zero but very small probability that firms extract resources at a loss, due to a lack of information 
or miscalculations of exploitation costs. 
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Figure 10: Left Progression of the extraction of units of a resource for a linearly increasing price Pit), where vertical slices of the distribution 
are gradually removed. Right Resulting supply of units as a function of time. 

strongly depends on the path taken by P(t). The magnitude of the flow of energy extracted depends directly on the rate 
of change of P(i). 

As an example, figure 10, left panel, depicts the progression of exploitation in time using a linearly increasing cost 
P(t). The black curve corresponds to the initial distribution, no. Each colour curve represents the number of units 
as extraction progresses after a certain time. These are equally spaced in time, which results from using a linear cost 
progression. The right panel of the figure shows the resulting supply of units as a function of time, calculated by taking 
the integral of the remaining number of units n(C,P(t),t) at each time value. 

However, when requiring a certain supply of units, the change in P(t) depends on the magnitude of no(C) at C values 
near that of P(t). For a large numbers of units situated narrowly in cost values C, the value of the market cost of 
extraction P(t) may be almost stationary. For low amounts of cheap units, the rate of change of P(t) will be very large. 
In this theory, therefore, the value of P(t) never endogenously decreases unless demand decreases. This is due to the fact 
that we have made abstraction of the effect of hoarding, and we thus assume that only the amount of energy intended 
to be immediately consumed is generated by producer firms. 

This also demonstrates the irreversibility in the process; if the price remains low, only low cost resources are used, 
while if the price increases, for example faster than consumption can occur, then significant amounts of low cost resources 
may remain in place while more expensive resources are consumed, and these two price behaviour assumptions could be 
chosen such that they generate exactly the same flow F. This corresponds to hysteresis and path dependence, akin to 
the change in entropy that occurs in irreversible processes in physical systems. 

S.1.7. The role of the parameter i>q 

The parameter i/q controls equation 6 by determining the fraction of reserves per unit time the system is allowed to 
consume during a unit of time (e.g. a year), for a multitude of internal reasons that need not be known in detail. All 
that is required to be known is that vq is indeed a constant of time, a fact fairly well demonstrated for oil and gas in the 
main paper. But what is the effect of changing the value of vq! This can be illustrated by the use of limiting values. 

For a value of i>q = 1, the complete amount of reserves can be consumed during a unit of time without requiring the 
price to increase. This results in the resources being exploited in perfect order of cost, with perfectly vertical slices of 
the distribution (as in figure 10) being consumed in order of cost. 

For a low value of vq, the opposite behaviour occurs, where a very small fraction of reserves at each value of C below 
P, producing a very low supply. This thus requires the price P to increase to values high enough that the supply meets 
the demand. At the extreme situation where the flow does not meet the demand, the price diverges to infinity and all 
resources of the distribution are exploited at equal rates. The slices of figure 10 thus become more or less horizontal. 

Thus it can be inferred intuitively that the role of the parameter vq controls the size of the reserves and how high the 
price is required to remain in order to supply the demand. It represents to some extent the rate at which resources can 
physically be taken out of the ground, but also to a certain degree how willing resource rich land owners are willing to 
exploit their resources instead of keeping them for a future where a higher price is expected. 
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S.1.8. Renewable resources 

This theoretical framework could also be used for renewable resources. In this case, n(C, t) would correspond to a 
distribution of resource producing units, and the price P(t) would be the price of electricity. Such a model would have 
quite different properties, stemming from the fundamental difference in the definition of n(C,t): it concerns units of 
flows of resources rather than resource units. Therefore, for a constant supply of resource, i.e. a constant flow, the 
level of resource use does not need to change, and thus the price P does not need to increase but simply to converge 
towards a constant value. This therefore shows that this model used for renewables would be very close to equivalent to 
a cost-supply curve framework, where a supply corresponds closely to a single cost value. There would be little gain in 
attempting to define such a model for renewables, and thus this is not done in the FTT model. 
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